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SUMMARY 

This report should be used as a reference guide when using the finite 
element method to approximate the static and dynamic behavior of flexible, 
rotating blades. Important parameters such as twist, sweep, camber, co-planar 
shell elements, centrifugal loads, and Inertia properties are studied. Compar- 
isons are made between NASTRAN elements through published benchmark tests. 

The main purpose of this report Is to summarize blade modeling strategies and 
to document capabilities and limitations (for flexible, rotating blades) of 
various NASTRAN elements. 


INTRODUCTION 

Rotating blades constitute vital components In several practical machines 
such as, turbofan and turboprop engines and helicopter rotors. They can fea- 
ture variable cross sections and spanwlse variations of sweep and pretwist. 
Because of their complex geometry, the finite element technique Is a natural 
choice for structural dynamic analysis. 

Until the mld-1960's, most of the structural dynamic analyses were based 
on continuum beam methods such as Raylelgh-Rltz, Myklestad, extended Holzer 
and Galerkln (see sample works of Rao and Carnegie (refs. 1 and 2)). In these 
methods the blade Is Idealized as a pretwisted cantilever beam. Hence, they 
are mainly useful for high aspect ratio blades. With Increased geometric com- 
plexities of modern blades, finite element methods have become more commonplace 
for static and dynamic analyses. One widely used general purpose finite ele- 
ment program Is NASTRAN. 

To employ finite element methods for rotating blades, analysts have to use 
extra care to obtain reliable results. Specifically, this Involves an Investi- 
gation Into element selection, grid modeling techniques, boundary conditions 
and effects of rotation. The purpose of this manual Is to discuss applications 
of the finite element technique to rotating blade structures. Typical test 
problems are examined and solution procedures are reviewed and presented. 
Emphasis Is directed towards NASTRAN. 


BACKGROUND 

The finite element technique offers many options for the structural 
dynamic analysis of swept and twisted blades. It Is common practice to use 
plate or shell elements of rectangular or triangular shapes. For vibration 
analysis of skewed cantilever plates, an early application of this technique 
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was used by Dawe (ref. 3) and Anderson, Irons and Zlenklewlcz (ref. 4), using 
rectangular and triangular elements, respectively. Olson and Llndberg (ref. 5) 
employed the same technique for finding free vibration characteristics of 
highly cambered blades using rectangular cylindrical shell elements. Rawtanl 
and Dokalnlsh (refs. 6 and 7) reported static bending and free vibration 
results of twisted cantilever plates. 

Presently, a joint government/industry/university research effort Is being 
conducted to compare finite element methods for twisted cantilevered plates to 
experimental data. The study Involves two phases: Phase I considers only non- 

rotating, twisted plates (refs. 10 to 12) and Phase II extends this effort for 
rotating plates (presently under Investigation). What has evolved from the 
first phase of this study Is that a wide variation In submitted finite element 
analyses exists and a careful application of this method Is needed to accu- 
rately predict the vibratory characteristics of highly twisted cantilever 
plates . 

Many of the finite element results reported In references 10 to 12 show 
mark disagreements In predicted natural frequencies, particularly for highly 
twisted plates. The triangular plate finite elements (Including COSMIC/NASTRAN 
CTRIA2 models) were typically deficient as were half of the shell and 
three-dimensional solid element formulations. Models using SAP quadrilateral 
plates and MSC/NASTRAN CQUAD4 shells, showed good convergence over a wide range 
of geometric parameters. 

Improved multiblade propeller engines, called propfans, are being devel- 
oped. Unlike conventional propellers, propfan blades are thin, flexible, shell- 
like structures with variable sweep and twist along their span. They are 
subjected to relatively high rotational speeds that can cause large deforma- 
tions at the blade tip. Because of the complex geometry of propfan blades, 
finite element methods are a natural choice for deformation analysis. However, 
recently documented finite element analyses of propfan blades have shown signif- 
icant disagreements with experimental data of modal frequencies (ref. 13). 

This report Is the first of two that deals with the finite element anal- 
ysis of flexible blade structures. A second report (refs. 44) Is a primer 
that summarizes the use of NASTRAN's large displacement and normal mode anal- 
yses for flexible, rotating blades. In the present report, previously pub- 
lished static and dynamic benchmark test problems are documented as they relate 
to blade configuration parameters such as twist, sweep and camber. The report 
Is organized as follows. A discussion Is first given of the types of approxi- 
mations likely to exist In a finite element analysis of practical blades. For 
completeness, applicable COSMIC and MSC/NASTRAN plate/shell and solid elements 
are briefly described. Additional topics related to the subject are examined, 
which Include element comparisons, mesh arrangement and refinement, co-planar 
shell elements, Inertia properties, and centrifugal forces. Numerical solu- 
tions and graphical comparisons, as reported by various Investigators, are given. 
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FINITE ELEMENT APPROXIMATIONS OF BLADE STRUCTURES IN NASTRAN 

The structural analysis of flexible, rotating blades has evolved from 
simple beam theories to sophisticated finite element techniques. Finite ele- 
ment grids for these blades can be divided Into two categories, two- 
dimensional and three dimensional. A two-dimensional finite element mesh links 
grid points on the blade's mid-surface by connecting these points with either 
flat plate or shell elements. A three-dimensional mesh Involves more than one 
surface of grid points, with two or more nodes In the thickness direction. The 
grid points are connected with "solid" or "brick" elements. 

A finite element analysis of rotating, arbitrarily-shaped blades using 
plates/shells or solids Introduces four different types of approximations, as 
quoted from Clough and Wilson of reference 14. 

"First Is the geometric approximation Involved In replacing the actual 
continuous structural surface by an assemblage of discrete structural elements" 
(ref. 14). For modeling of swept and twisted blade surfaces In NASTRAN, sev- 
eral structural elements are available In the form of flat triangles (TRIA2, 
TRIA3.TRIA6) , flat rectangles ( QUAD2 , QUAD4 , QUAD8) or hexahedrons (HEX8,HEX20, 
HEXA8.HEXA20) . Each of these elements possess their own special features that 
are pertinent to rotating blade structural analysis (refs. 8 and 9). 

"The second basic assumption Is that the strains (and/or displacements) In 
each element may be only of limited form, as specified by a set of nodal Inter- 
polation (or "shape") functions" (ref. 14). Depending on the actual stress 
variation placed on the blade at a given rotational speed, these assumed dis- 
placement orders may provide a reasonably good approximation. The NASTRAN 
higher order elements (TRIA6,QUAD8,HEXA20,HEX32) that have quadratic or cubic 
displacement variations, may provide better results than the lower order ele- 
ments for a given degree of mesh refinement. On the other hand, more refined 
elements that feature higher order displacement fields with additional global 
degrees of freedom (DOF) would not be suited for the format of NASTRAN's gen- 
eral element. A practical "rule of thumb" to consider In the analysis Is that 
"simplest can be best." That Is, theoretically, any Inherent error found In 
the calculated results of a coarse mesh consisting of simpler finite elements 
tends to vanish as the mesh Is refined. 

"A third approximation which may be present In some finite element anal- 
yses Is that the assumed displacement patterns may not maintain Interelement 
compatibility as the blade structure Is loaded and deformed" (ref. 14). It Is 
quite apparent that the problem of "co-planar elements" arises when using 
NASTRAN's TRIA and QUAD shell elements In a large displacement analysis 
(NASTRAN SOLUTION 64 — S0L64). This element breakdown comes from a variety 
of sources. For Instance, when modeling blades with high sweep and/or twist, 
the chord dimension at sections along the span are such that many elements may 
have to be used to achieve a satisfactory level of solution accuracy. Element 
Incompatibility can develop when neighboring elements meet at a finite angle 
during any particular Iteration of S0L64. This causes the analysis to abort at 
the global level, because the shell elements have 5 local degrees of freedom 
(DOF) per node and the global transformation equations utilize 6 global DOF 
per node. 

Finally, the fourth approximation Is that the three-dimensional blade sur- 
face may be treated as a two-dimensional surface with simplifying constraints 
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on displacement variations through the blade thickness. In COSMIC/NASTRAN ‘ s 
TRIA2 and QUAD2 element formulations, the Klrchoff thin plate theory of plane 
sections Is assumed. However, less restrictive assumptions of forming a shell 
element as a degenerate three-dimensional solid element can be advantageous. 
What results Is that shear distortions may be accounted for without difficulty 
by this approach. Ahmad (ref. 16) used this approach by Introducing a con- 
straint that lines through the element thickness displace In translation or 
rotation without distortion. Zlenklewlcz (ref. 15), Iron (ref. 16), and Pawsey 
(ref. 17) also used degenerate solid elements to analyze arbitrary thin shells. 
By using modified Integration points to calculate shear strain energy terms In 
the solid elements, overly-stlff element behavior In bending was suppressed. 

From a practical standpoint, a good finite element Is not necessarily the 
one which exhibits a monotonic convergence. A more Important criterion Is that 
the deviation from the exact result should be small, even when the blade Is 
subdivided In a relatively coarse mesh. This point Is especially Important 
when using more refined shell and solid element assemblages, since they require 
a great deal more computer memory and cost. 


DESCRIPTION OF NASTRAN'S PLATE/SHELL AND SOLID ELEMENTS 

From the preceding discussion, It Is apparent that the type of element 
utilized may have an effect on solution accuracy and rate of convergence. For 
completeness, this section presents, as a brief catalog, those NASTRAN plate/ 
shell and solid elements that may be used. The elements are shown In figures 1 
to 10. Their principal features are listed below. 

COSMIC/NASTRAN plate elements have two uncoupled stress systems, l.e., 
membrane and bending. COSMIC/NASTRAN solid elements are based on a standard 
Isoparametric formulation, where Identical approximation functions are used for 
the element geometry and displacement. Differential stiffness formulations are 
used for both the plate and solid elements In COSMIC/NASTRAN. These elements 
are described below: 

TRIA2 Is a constant thickness triangular element with both In-plane and 
bending stiffness (fig. 1). The stiffness formulation Is a superposition of a 
constant strain membrane triangle with linear displacement order and of a basic 
bending triangle with an eight term, cubic normal displacement order. Important 
details of the stiffness derivation Include: (1) bending slopes of the element 

are obtained from the definition of transverse shear strain (assumed constant 
In the element), and (2) a constraint such that the bending slope about the 
x-axIs varies linearly from side a-b (fig. 1). Two TRIA2 elements Joined along 
this boundary will have continuous displacements and slopes at all points. 

Both lumped and consistent mass formulations can be used. 

QUAD2 Is a quadrilateral element comprised of two sets of overlapping tri- 
angular elements (fig. 2). The a-b sides of the basic NASTRAN triangle lies 
along a diagonal to ensure displacement and slope continuity within the ele- 
ment. Since membrane and bending stiffness coupling Is not available In 
COSMIC/NASTRAN, the QUAD2 Is formed by adding the effects of the overlapping 
quadrilateral membrane element (QDMEM) with those of the overlapping quadri- 
lateral bending element (QDPLT) . 



OTHER PLATE/SHELL ELEMENTS Alternative element formulations based on pre- 
vious analyses posted In the literature can be used In COSMIC/NASTRAN. For 
example, one might consider using the Clough bending triangle (ref. 18), which 
can be superimposed with a NASTRAN membrane triangle to form a flat shell ele- 
ment (fig. 3). However, this element Is archaic for thin, flexible blade 
applications In that It Imposes additional constraints to ensure displacement 
continuity within the element. Consistent mass formulations of this element 
produce vibration frequencies on the high side, because constraint approxi- 
mations, used In the stiffness calculation, are progressively applied In the 
mass formulation. 

A second choice might be COSMIC/NASTRAN ' s TRSHL element. This element Is 
a higher order triangular shallow shell element. The element approximation 
comprises a quadratic membrane and a qulntlc normal displacement function. 

Usage of TRSHL Is low, most likely as a result of Its Inferior convergence on 
shell benchmark test problems. 

There Is a more attractive alternative to using QUAD2 elements In COSMIC/ 
NASTRAN. Because of Its potential for modeling curved surfaces, a quadri- 
lateral shell element can be formed (manually In NASTRAN) out of four subdomain 
triangles with corner nodes not In the same plane (see fig. 4). The a-b side 
of the basic NASTRAN triangle Is placed on an exterior edge to ensure Interele- 
ment continuity. The Interior point c Is located at the Intersection of 
straight lines connecting the midpoints of the sides (or at the average of the 
four corner node coordinates). The membrane and bending DOF at point c can be 
condensed out of the stiffness equations. Thus, the quadrilateral effectively 
has only 20 DOF per element. Although static condensation of the Interior 
node Introduces an additional constraint on this element. Its level of engi- 
neering accuracy Is quite good. Similar elements proposed by Clough and 
Fellppa (ref. 19) and Clough and Johnson (ref. 20) were placed In the SAP gen- 
eral purpose finite element program. The twisted plate results shown In 
reference 10 using these elements were shown to perform well against those of 
other finite element formulations. 

HEX8 and HEX20 are standard eight and twenty noded Isoparametric brick 
elements (figs. 5 and 6). The element geometry and displacements are approxi- 
mated by linear and quadratic functions, respectively In the three local direc- 
tions. The element coordinate systems are defined on the CHEXA Bulk Data 
cards. Material reference Is specified on the PHEX card. Reduced order 
Integration Is available (via NGP-fleld on PHEX Bulk Data card) for both HEX8 
and HEX20 . 

The MSC/NASTRAN element library Includes both flat and curved shell ele- 
ments, which support both membrane and bending action. The analyst has a 
choice of four control options (via MID fields on PSHELL card) for stiffness 
formulation: (1) membrane-only behavior, (2) bending-only behavior, (3) com- 

bined bending and membrane behavior, and (4) materially coupled membrane-bending 
behavior. These stiffness control parameters are Important options for struc- 
tural analyses. The shell elements In MSC/NASTRAN Include the following: 

TRIA3 and QUAD4 are the three noded triangular and the four noded quadri- 
lateral elements, shown In figure 7 and 8, respectively. General capabilities 
Include: (1) Identical approximating functions for element geometry and dis- 

placements, (2) 5 DOF per node, l.e., two membrane displacements and one bend- 
ing displacement and two bending rotations, (3) thickness can vary over the 
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element surface, (4) lumped and consistent mass formulation, (5) selective 
Integration for shear terms (QUAD4 only), (6) enforcement of bending curvature 
compatibility (QUAD4 only), (7) transverse shear flexibility to account for 
bending behavior (QUAD4 only), and (8) full geometric nonlinearity. All forces 
and moments are evaluated at the centroid of the element. Use of QUAD4 ele- 
ments requires that all Interior angles be less than 180°. A special feature 
of the QUAD4 elements Is that all points do not have to lie In one plane. When 
modeling highly swept and twisted blades, enough QUA04 elements should be used 
so that Individual element warping Is small. According to (ref. 9), "for 
modeling curved shell structures a tentative rule Is that the Included angle 
per element should be 10° or less to obtain errors of 4 percent or less In 
deflections. 

MSC/NASTRAN ' s curved shell elements are the following: 

TRIA6 and QUADS are the six noded triangular and eight noded quadrilateral 
curved thin shell elements, respectively shown In figure 9 and 10. The assumed 
displacement and geometry of these elements are of higher-order than their low- 
order cousins, described above. Also, these curved shell elements can connect 
six and eight grid points, respectively and they contain all of the features 
of their flat shell counterparts. Note that QUADS uses a standard Isoparam- 
etric formulation (l.e.. Identical approximation functions are used for both 
element geometry and displacements). QUAD8 does not have transverse shear 
flexibility terms like Its QUAD4 counterpart. For the TRIA6 and QUAD8 elements, 
all forces and moments are evaluated at the centroid and at the vertices of the’ 
element. 

MSC/NASTRAN 1 s solid elements are the following: 

HEXA8 and HEXA20 The HEXA series of elements (figs. 5 and 6) In MSC/NASTRAN 
are Improvements to the HEX8 and HEX20 six-sided brick elements In COSMIC/ 
NASTRAN. These elements use Isoparametric representation to characterize the 
element behavior. The HEXA8 Is the three-dimensional analogy of the two- 
dimensional bilinear quadratic QUAD4 element and the HEX20 Is the three- 
dimensional extension of the two-dimensional biquadratic cubic QUADS element. 

The HEX8 element contains the same type of shear correction used to Improve 
the membrane shear behavior of Its two-dimensional cousin. In addition, a 
reduced order of Integration (via IN-fleld on PSOLID card) can be used for the 
shear strain contributions to the stiffness calculations. This Is discussed 
In the next section. 


STATIC AND DYNAMIC BENCHMARK TESTS FOR BLADE STRUCTURES 

The proceeding element descriptions are Insufficient to make a decision 
on which elements are best to use for a particular blade. Lacking In the 
NASTRAN documentation are the performance of Its elements In standard benchmark 
test that are applicable to swept, twisted or cambered blades. For example, 
some test problems that measure efficiency of plate/shell or solid elements’ 
Include: (1) Scordells-Lo shell problem, (2) patch test (which measures con- 

vergence of an assemblage of elements), (3) untwisted cantilever beam with 
point load at the free end, (4) swept cantilever beam with point load at the 
free end, (5) swept plate with uniform load, (6) twisted cantilever beam with 
point load at the free end, (7) twisted cantilever plate (normal mode vibra- 
tions), and (8) cantilever blade with camber (normal mode vibrations). Many of 



these tests have been documented In surveys of finite element codes, namely 
Fong and Jones (ref. 21) and Harder (ref. 22). 

Although the results of these published tests do not completely test the 
elements, they are sufficient for making preliminary decisions regarding finite 
element modeling of practical blades. More Importantly, they reveal Inherent 
deficiencies In a number of the NASTRAN plate/shell and solid element formula- 
tions. The succeeding subsections examine these benchmark test problems. 
Specific blade modeling aspects are discussed such as Irregularity of mesh 
arrangement, mesh refinement, element length and thickness ratios, element 
warp and skewness, and modified Integration within elements. 


The Scordells-Lo Shell 

Figure 11 shows the problem and table I and figure 12 display the conver- 
gence behavior of NASTRAN QUAD and HEXA elements In a rather severe test of 
bending behavior. The data has been compiled from references 21 and 22. 

Figure 12 shows excellent convergence of QUAD4 between the analytical shallow 
shell and the exact deep shell solutions. This suggests that this element Is 
quite suitable for highly cambered blade models. Although the results of TRIA2 
tend to overshoot the mark, the convergence Is Inconclusive. With a coarse 
mesh this Is to be expected, since the transverse shear strain Is assumed con- 
stant In the TRIA2 element formulation. QUAD2 are typically overstiff even 
for large numbers of DOF. Transverse shear flexibility In QUAD4 explains Its 
obvious suitability for this benchmark test problem. 

To shed light on NASTRAN's higher order elements, consider the normalized 
vertical deflection data for the Scordells-Lo shell In table I. This data 
reveals that QUAD4 and HEX8 exhibit a strong convergence from above, and that 
QUAD8 and reduced Integration forms of HEX20 show an oscillatory pattern of 
convergence. Clearly, QUADS portrays the best numerical behavior for this 
problem with Its quadratic shape functions for Its geometry and displacements. 
QUAD2 and HEX20 (without reduced Integration) are typically overstiff. 


The Patch Test 

In some displacement-based element theories approximating functions may 
violate Interelement continuity of displacements. This causes Infinite strains 
along the element boundaries. Elements characterized by this behavior are 
called "nonconforming elements." In the limit the formulation may converge to 
the correct answer with Increasing mesh refinement. Necessary conditions are 
(1) that zero strain exists In the element during rigid body motion and (2) 
that the assumed displacement function contains terms for which constant strain 
modes are obtained. 

To ensure that convergence may be achieved, nonconforming elements must 
pass a patch test (ref. 40). This test Involves an arbitrary assemblage of 
oddly shaped elements subjected to a constant strain/stress condition. If the 
resulting nodal displacements cause a constant stress state within the element, 
then the patch test Is not violated. Nonconforming elements that pass this 
test have convergent results with mesh refinement, and at times will exhibit 
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better behavior than conforming elements. However, note that nonconforming 
elements that violate the patch test will not necessarily produce unacceptable 
results . 


Figures 13 and 14 depict typical patches for plate and solid elements In 
a severe test for bending action. Results tabulated In table II (ref. 22) 
show maximum error In stresses for QUAD and HEXA elements. Some observations 
can be made from this data. HEX8 with linear strain and HEX20 with quadratic 
strain are conforming elements. As expected, the patch test Is not violated. 
However, the HEX20 with reduced Integration, which Is nonconforming, passes 
the test. 

Of the plate elements, QUAD2 falls the patch test, because bending rota- 
tions In this element are derived from a constant transverse shear strain var- 
iation within the element. Thus, It cannot represent a constant curvature 
state. On the other hand QUAD4 exhibits both force and moment convergence, 
mainly due to correction terms for transverse shear flexibility Imbedded In 
the formulation. These correction terms do not change the rigid body modes of 
the element nor do they eliminate the pure shear deformation behavior of the 
element. They only modify the magnitude of the stiffness coefficients and the 
locations where the strains are computed. QUADS Is an Isoparametric, conform- 
ing element without the above devices to Improve Its bending behavior. This 
may account for Its poor convergence In table II. 


Effect of Element Distortion 

Careful attention should be given to preserving Isosceles, square and cube 
shapes of NASTRAN TRIA, QUAD and HEX elements, respectively. To Illustrate the 
effect of distortion In NASTRAN elements, figure 15 displays three untwisted 
cantilever beams with different mesh arrangements. The geometric dimensions, 
material parameters and load conditions are shown In the figure. In the first 
case the elements are undistorted. This case Is designated as a regular mesh 
(Reg) arrangement. For the second case the elements sides are skewed at a 
45° angle to form trapezoidal shape elements In the first of two Irregular mesh 
(IRReg) arrangements. Finally, the last case shows an IRReg arrangement of 
parallelogram shape elements with sides also skewed at 45°. 

Normalized displacement response at the beam tip Is tabulated In table III 
for extenslonal , flapwlse bending, edgewise bending and torsional modes. Some 
observations can be drawn from the numerical data. First, the extenslonal and 
torsional modes are Independent of element distortion. Identical results were 
obtained for the regularly and Irregularly-shaped elements. This Is Indicated 
In table III by the (BOTH) notation. Second, the change In results between the 
two IRReg arrangements Is negligible. Third, MSC's QUAD4 seems to show a mark 
degradation of accuracy In the flapwlse bending mode of IRReg models, while Its 
three-dimensional cousin, HEXA8, exhibits an error In both the flapwlse and 
edgewise bending modes. In a paper on QUAD4 ' s formulation (ref. 23), the 
author cites this effect of element distortion: "After the paper was submitted 

for publication, It was discovered that large errors occur when the skew angle 
of the element exceeds 20°. This error was traced to coupling between trans- 
verse shear strains, and has subsequently been corrected" (1978). Apparently, 
this deficiency Is still present In QUAD4 as Illustrated by these more updated 
results (ref. 22), and thus, the 20° limit on element skew should be used when 
using QUAD4 or HEX8. 
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With the exception of the above remarks, all of the elements In table III 
display good to excellent convergence. Note that QUAD2, although not shown, 
performed poorly In the edgewise bending and torsional modes (ref. 22). This 
deficiency Is, again, related to the constant transverse shear strain assump- 
tion In the QUAD2 element formulation. 


Cantilever Beam with Sweep 

Recently developed propfan blades are highly swept at the tip (approxi- 
mately 35° In some proposed models (ref. 33)). To Illustrate NASTRAN plate/ 
shell and solid element performance on swept blade configurations, consider 
the curved cantilevered beam shown In figure 16. The beam Is curved (or swept) 
In Its own plane with an Included angle of 90°. Table IV displays results for 
two tip load conditions. As expected, the QUA02 exhibits overstiff results for 
both directions. But, QUA04 and HEX20 both show average convergence In the 
In-plane response direction and good convergence In the out-of-plane direction. 
Better still, reduced Integration forms of HEX20 and QUAD8 portray the best 
results by far. 


Morley's Simply-Supported Swept Plate (ref. 24). 

In June 1983, a user project on sweep effects In plates was proposed In 
Finite Element News (ref. 25). The Investigators received participation from 
numerous researchers and users. The scope of submitted finite element results 
revealed 33 plate bending element formulations and employed 19 finite element 
computer codes. Three finite element mesh densities were requested to test 
convergence of predicted center point displacements of simply-supported plates. 
A total of four sweep angles were requested to show distortion sensitivity of 
the finite elements. This section discusses the MSC/NASTRAN results. 

Figure 17 shows a schematic description of the proposed plate model. The 
series solution developed by Morley (ref. 24) Is given In table V. The MSC/ 
NASTRAN results reported (ref. 25) are for the QUAD4 and TRIA3 elements. The 
triangles are chosen so that all rectangles are cut along the shorter diagonal, 
as shown In figure 17. At zero sweep, the plate mesh Is symmetrical. This 
ensures a maximum number of triangles at the center point of the plate, where 
maximum force and response Is anticipated. 

Numerical data tabulated In tables VI and VII show the percentage error 
In calculated moment and response of using QUAD4 and TRIA3 elements. These 
percentage errors are based on Morley's series solution for the maximum and 
minimum bending moments and the vertical displacement at the plate center. 
Results for the three mesh conditions are given for sweep angles of 30 to 90° 
as designated In figure 17. Although the tabulated results show the effects 
of sweep, they also Include the effect of aspect ratio (a/h) (see fig. 17). 
Since the plate has a constant boundary length (a), the value of h will change 
for each sweep angle. The aspect ratio for an unswept plate Is assumed to be 
unity. The aspect ratio values Increase with sweep angle to a final value of 
two for 30° sweep. 
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Additional points can be drawn from the data In tables VI and VII. Gen- 
erally, with both QUAD4 and TRIA3 element assemblages, an error Is apparent 
for plates with sweep angles between 40 and 60°. As sweep angle Is Increased, 
the difference In accuracy of calculated moment and response between the QUAD4 
and TRIA3 elements diminish. Intuitively, this Is suspected to be an effect of 
aspect ratio (which Is further discussed In a later section). For the larger 
sweep angles at a given mesh densities, the QUAD4 and TRIA3 elements become 
quite distorted In shape. This distortion Is a major contribution to the error 
of the predicted results. Thus, It Is safe to suggest that a refined mesh 
yields more accurate results for highly swept plate models. This will yield 
better-shaped QUAD4 and TRIA3 elements In the mesh. 


Cantilever Beam with Pretwist 

A cantilever beam with 90° twist from root to tip Is shown In figure 18. 
Table VIII compares NASTRAN triangle and solid element prediction of displace- 
ment response at the beam tip for two bending directions. All elements per- 
formed quite well for each response mode with TRIA6 and reduced Integration 
form of HEX20 as particular standouts. 


Element Aspect Ratio and Modified Integration 

An Indirect effect of element aspect ratio Is shown by the normalized 
deflection response data for a quarter-plate mesh of NASTRAN's QUAD or HEX 
elements In figure 19 and table IX. Plate aspect ratios of 1 and 5 are used 
as example test cases. Note that the HEX20, shows poor results for the very 
coarse meshes. However, the HEX20R Improves these results very significantly. 
The HEX20R stiffness coefficients are calculated by modifying the evaluation of 
strain energy In the HEX20 . The essential difference here Is that the number 
and location of the Integration points used within the element to evaluate the 
shear strain energy terms are chosen differently. Details of this technique 
are Inappropriate here, however, the basic Idea can be explained by referring 
to the Isoparametric quadrilateral In figure 20. 

It can be seen (fig. 20(b)) that an actual residual bending occurs due to 
transverse shear action, when the element Is subjected to a pure bending 
stress. This residual bending becomes very significant as the aspect ratio of 
the element Is Increased. The assumed displacement functions require that the 
element edges remain straight under pure bending (solid lines), where In fact, 
the edges should be curved (dashed lines). If the transverse shear effect were 
omitted the element stiffness representation would be substantially Improved. 
Note that the strain due to pure shear stress (fig. 20(c)) should remain. A 
practical procedure to obtain the shear strain terms for HEX elements Is to 
use Its centroid as the Integration point; this neglects the shear strain of 
figure 20(b), yet retains that of figure 20(c). 

The use of modified Integration points, chosen so as to suppress undesir- 
able element behavior, Is noteworthy. Particularly, It Is quite effective In 
beefing up the bending action of three-dimensional Isoparametric elements, when 
they are used to model thin or moderately thick blade surfaces. This element 
modification Is especially beneficial as the element aspect ratios becomes 
large. Note that this element Improvement device Is not fall-safe. Because, 
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as element aspect ratios become too large, mesh redefinition may be the only 
alternative approach. 


Cambered Blade and Twisted Blade Vibrations 

This section examines benchmark tests for vibration analyses. The cam- 
bered compressor blade described In figure 21 was analyzed by MacNeal (ref. 23) 
In a severe vibration analysis test of QUAD4 elements. The numerical results 
are tabulated In table X. Eight free vibration modes are given along with some 
comparative data from Olson and Llndberg (refs. 5 and 29). All modes were pre- 
dicted extremely well. However, the second mode Is a particular stand-out, 
because Its refined results agrees more with Olson's analytical calculations 
and not with experiment. The probable source of error comes from support com- 
pliance In the prediction of the second mode. Because this mode Is typically 
characterized by no surface nodal lines In camber blade configurations. It Is 
difficult to accurately model Its support condition. 

Although finite element analysis for plate vibrations has become common- 
place, "significant differences In the published results for various methods 
of analysis have raised some doubt concerning the adequacy of these methods to 
accurately predict the vibratory characteristics of highly twisted cantilever 
plates," according to Klelb, et al. (ref. 10). They continue to explain, "An 
example of the types of difference existing In the literature Is the predicted 
first bending frequency for pretwisted plates. The predicted dependency of 
frequency on twist angle for twisted, cantilever plates spans the spectrum from 
Increasing significantly to decreasing significantly. To be more specific, 
the trends seen In the literature are: reference 30 - significant Increase, 

references 31 to 33 - slight Increase, reference 34 - no change, references 35 
and 36 - significant decrease, references 37 and 38 - slight decrease and then 
an Increase, and references 7 and 39 - significant decrease." 

Two proposed TRIA2 and QUAD4 grids for the twisted plate model In 
figure 22 were recently submitted to a finite element users project (refs. 10 
to 12). The TRIA2 grid Is shown for a plate of aspect ratio (a/b) equal to 1 
and a total plate twist of 45°. A total of 360 DOF and 128 elements are used 
with 9 nodes In the spanwlse and chordwlse directions. A lumped mass formula- 
tion Is utilized. The QUAD4 grid represents a plate with a/b = 3 and a total 
twist of 60°. A total of 550 DOF and 100 elements are used with 11 nodes In 
the spanwlse and chordwlse directions. A consistent mass formulation Is 
employed. Figure 22 also shows the dependence of frequency parameter on twist 
angle for a/b = 1 only. The results of the TRIA2 and QUAD4 elements and those 
of experiment are highlighted over the finite elements results used by other 
Investigators. It Is apparent that the predicted first bending frequency var- 
iation Is significantly Increasing for the TRIA2 model (highlighted with dots), 
whereas the trend for the QUAD4 model Is significantly decreasing (highlighted 
with stars). Note that the two experimental results (shown as dotted lines) 
depict a significantly decreasing trend. In all modes shown, the TRIA2 models 
Is predicting overly-stlff results compared to those of the QUAD4 model. 
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Summary of Test Results 


An attempt has been made to grade the TRIA, QUAD and HEX elements In their 
performance In the above tests. The results are given In tables XI and XII. 

An average rating and an approximate grade for each element Is compiled at the 
bottom of the table. Based on the judgments of the author, QUAD4 and TRIA3 
exhibited the best overall performance and Is the recommended choice amongst 
the NASTRAN plate elements for practical blade applications. Both of these 
elements pass the patch test, which ensures (In the limit) convergence of solu- 
tion. Of the NASTRAN solids examined, HEX20 with modified Integration per- 
formed the best and Is highly recommended for three-dimensional analysis of 
thin or moderately thick blade planforms. 


CO-PLANAR SHELL ELEMENTS IN NASTRAN 

While using NASTRAN shell elements, undesirable difficulties arise when 
these elements are joined at finite angles. This element joining Is commonly 
referred to as "co-planar elements." The problem Is especially acute, while 
the blade Is loaded and deformed during particular subcases of NASTRAN SOLUTION 
64 (large displacement analysis). Essentially, a surplus of degrees of freedom 
(DOF) corresponding to an In-plane rotation originates In the transformation 
phase of NASTRAN shell elements and can cause singularity of the global stiff- 
ness matrix. What results Is a solution that either requires a large number 
of subcases for convergence or leads to nonconvergence. 

To clarify this point, consider the example shown In figure 23. Here, 
representative triangular elements have Individual nodal parameters correspond- 
ing to membrane and bending actions. If the triangles are not co-planar In an 
assembly, the common node A Is capable of resisting 6 Independent responses, 
when the nodal parameters are transformed Into global coordinates. However,' If 
the elements are co-planar, only 5 DOF will exist and thus, global transfor- 
mations will create zero stiffness coefficients In the sixth DOF. The TRIA3 
and QUAD4 elements are more likely to have co-planar difficulties than the 
TRIA6 and QUADS, since the midside nodes ensure better compatibility of the 
local normals of adjacent elements. 

To overcome this source of difficulty, a variety of practical procedures 
have been used In the analysis of arbitrary shells. An Indirect approach 
Involves taking, as the normal, the average of the normal directions associated 
with each of the elements connected at a particular node of the shell. When 
the element stiffness has been transformed to the global system, the sixth DOF 
Is eliminated by deleting the corresponding row and column from the global 
stiffness matrix. Eliminating this DOF Is equivalent to constraining the shell 
structure against this motion. 

Physically, the shell Is quite stiff In this normal rotation, so the 
Imposed constraint has a negligible effect on Its response. This point Is 
clarified by the data In table XIV. A cylindrical shell Is clamped at one end 
(fig. 24), and allowed to deflect under Its own weight (ref. 41). Five differ- 
ent lengths are considered, ranging from 25 to 200 ft. Plane quadrilateral 
shell elements are used with meshes having 5 nodes In the circumferential 
direction and 6 to 11 nodes In the spanwlse direction, depending on the shell 
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length. The computed vertical deflections at the free end show that for the 
very flexible shell configuration, the effect of the sixth DOF constraint Is 
negligible. 

A more direct remedy to the co-planar element problem Is to Insert an 
arbitrary stiffness coefficient In the sixth DOF of the element stiffness 
matrix. After transformation to the global system, a perfectly well-behaved 
set of equations Is achieved from which, all displacements and stress result- 
ants, Including those corresponding to the sixth DOF, are obtained. A ficti- 
tious set of normal rotation stiffness coefficients for triangular shell 
elements (table XV) have been suggested by Zlenklewlcz (ref. 42). Similar 
matrix terms can be added for rectangular elements. These additional terms are 
far better than just placing small terms on the diagonal, because they pass the 
patch test, when added to the element stiffness matrix. Like the first method, 
the effect of these fictitious springs on calculated results Is negligible. 

A new parameter card K6R0T exists In MSC/NASTRAN (versions 63 and higher) 
that does this automatically. Actually, for QUAD4 and TRIA3 only, the added 
artificial stiffness coefficients are similar to those of reference 42. A 
recently completed study (ref. 43) revealed that values of the parameter K6R0T 
of 1000, 10 000 and 100 000 alleviated the singularity problems and gave 
smoother responses for advanced turboprop blade models on an element-toelement 
and node-to-node basis. 

Figure 25 displays strain gage locations and two finite element grids 
(using QUAD4 and TRIA3) for a SR-3C-3 propfan blade, which was analyzed In the 
study of reference 43. The value of K6R0T was varied to show that frequencies 
and response are not significantly changed with the variation of this arti- 
ficial plate normal stiffness. It was also noted that predicted strains were 
not significantly affected by changes In K6R0T. Additionally shown In the 
study was that the element-to-element strain variations became much "smoother" 
when the triangles were made more nearly equilateral or when the triangles 
were eliminated and QUAD4 elements were used (fig. 25). 

Table XVI presents frequency results for the QUAD4 grid model at 0 rpm 
using K6R0T = 0, 10 000 and 100 000. A faster convergence rate to a steady- 
state response of the blade Is observed with K6ROT = 10,000, where a total of 6 
subcases were needed In NASTRAN Solution 64, Instead of 25 subcases (with 
K6R0T = 0). Additionally, Increased values of the K6R0T parameter resulted In 
small changes In the calculated modes. NASTRAN runs using K6R0T = 10,000 with 
TRIA3 are also shown In table XVI. Strain gage results for the propfan blade 
rotating at a speed of 8508 rpm are given In table XVII. Strain gage locations 
are shown In figure 25. Gage 1 measured a strain In the radial direction of 
the blade. At gage 2 a strain was measured In the circumferential direction, 
while a shear strain was measured at gage 3. Again, measured strains at the 
gage points of both the QUAD4 and TRIA3 elements were not significantly 
affected by changes In the K6R0T parameter. 

From the study (ref. 43) It was concluded that the use of K6R0T to add 
artificial plate normal stiffness significantly reduced element-to-element 
strain variations In the finite element models of the flexible blades examined. 
Additionally the nonlinear steady state solution converged much faster than 
when alternative procedures were used. A value of K6R0T = 10,000 was shown to 
give good results for the practical blade models. 
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A companion memorandum (ref. 44) gives specific NASTRAN data Instructions 
on the use of PARAMeter card K6R0T. 


EFFECT OF CENTRIFUGAL STIFFENING AND SOFTENING FORCES 

Consider a straight cantilever blade undergoing bending vibrations per- 
pendicular to the plane of rotation. The centrifugal tensile force tends to 
stiffen the transverse bending elastic springs. In vibration analysis, this 
force generally Increases the square of the natural frequencies In proportion 
to the square of the rotation speed. The constant of proportionality Is known 
as the Southwell coefficient (ref. 1). Centrifugal stiffening effects are 
placed In the finite element formulation through supplemental strain energy 
terms due to the Initial stress of steady-state rotation. 

In a rotating reference frame, there Is an additional force that acts In 
the radial and circumferential directions of the blade. As a mass point dis- 
places radially outward, there Is an effective Increase In the total magnitude 
of centrifugal force carried by this point that Is proportional to Its total 
displacement. When this effective centrifugal force Increase Is algebraically 
added to the Incremental elastic stiffness (l.e., brought from the right-hand 
to the left-hand side of the equation of motion), a reduction In elastic stiff- 
ness of the system Is apparent. This stiffness reduction Is called “centri- 
fugal softening." In vibration analysis, this force has the effect of reducing 
the natural frequencies. Centrifugal softening effects are Incorporated In the 
finite element formulation by subtracting the mass times the square of the 
rotational speed from the stiffness matrix (see ref. 44 for detailed 
formulations) . 

To Illustrate the effect of the centrifugal stiffening and softening 
forces, consider the flat cantilever blade with length ratio equal to 3. The 
axis of rotation Is placed at the root section. This problem Is the second of 
two phases of a joint research effort on vibrations of twisted blades (refs. 10 
to 12). Tables XVIII and XIX each show five sets of submitted results of non- 
dimensional frequencies for the first five modes at three-nondlmenslonal 
rotation speeds (0,2,6). For the results In table XVIII, all mass points of 
the plate are spinning In the plane of rotation, whereas a 90° setting angle 
of the plate Is assumed for the cases of table XIX. 

As expected, an apparent Increase In the frequency parameter Is observed 
as a result of the centrifugal stiffening effect due to rotation. Investiga- 
tors A,C and D Included centrifugal softening In their formulations. Investi- 
gators B and E omitted this effect. From the numerical trends In table XVIII, 
It Is quite obvious that the centrifugal softening effect Is essential to the 
accuracy of frequency calculations of rotating blades. With a 90° angular 
rotation of the plate, the centrifugal softening effect disappears. Thus, this 
effect Is dependent on the blade's setting angle at the axis of rotation (or 
root section) . 

The question of Inserting centrifugal softening terms In a NASTRAN run of 
a practical blade Is Important, especially when large displacement effects are 
considered. NASTRAN specific data Instructions for Insertion of these terms 
are Inappropriate to discuss herein. These aspects are dealt with In 
reference 44. 
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EFFECT OF CORIOLIS FORCES 


In addition to the centrifugal stiffening and softening forces present In 
rotating blades, there exist Coriolis forces that affect the blade response at 
relatively high rotation speeds. Coriolis forces are produced when a mass 
point of the blade displaces relative to a rotating reference frame. These 
forces are perpendicular to the plane containing the rotation axis of the blade 
and the Instantaneous velocity vector of the mass point. They are proportional 
to the vibratory velocities, but unlike viscous damping forces, are path- 
independent In nature. In a finite element formulation, the contrast between 
Coriolis and damping forces Is the skew- symmetric versus symmetric form of 
their matrices. 

Many Investigators have considered Coriolis forces In their blade vibra- 
tion analyses and have commented on the significance and Importance of this 
effect. One such work Is that of Sreenl vasamurthy and Ramamurtl (ref. 47) In 
the effect of Coriolis forces on rotating plates. Here, triangular plate 
finite elements with uncoupled membrane and bending stress behaviors were 
utilized. Their results showed that for a uniform thick plate of aspect ratio 
equal to one and rotating at a speed near the first nonrotating frequency, the 
first bending and torsional frequencies decreased by 3.23 and 8.02 percent 
respectively, due to Coriolis forces. The overall effect of these forces were 
most significant for plate setting angles of 45°. The authors stated that for 
plates rotating at high speeds Inclusion of Coriolis forces In the computation 
of natural frequencies Is desirable. 

However, upper bound solutions for the free vibrations of rotating canti- 
lever plates have recently been reported by Co (ref. 48). For thin plates of 
aspect ratio equal to one, Co found the most significant change In the funda- 
mental mode due to Coriolis terms to be less than 0.1 percent at a rotation 
speed equivalent to three times the nonrotating frequency. It was also noted 
that by Increasing the plate thickness the Coriolis effect Increases, but no 
more than 1 percent for the fundamental mode. 

The Rltz method utilized In reference 48 provides upper bound solutions. 
Generally, the significance of Coriolis forces arises at very high rotation 
speeds of thicker blades mounted at high root setting angles. Modern propfan 
blades operate at more realistic rotation speeds and do not fall Into this 
category. Thus, In these blades Coriolis effects are typically small and 
negligible. 


LUMPED VERSUS CONSISTENT MASS MATRICES 

The mass matrix formulation In NASTRAN Is lumped unless the analyst 
requests the consistent formulation by means of a PARAM COUPMASS Bulk Data 
card . 


In a lumped mass formulation Inertia properties are assigned to the trans- 
lational and rotational DOF only. The total mass In each element Is distri- 
buted to the nodal points In an averaging-type of fashion, depending on the 
element type. In consistent mass formulations the mass matrix Is calculated 
with the same approximating functions used In the stiffness formulation. That 
Is why the word "consistent" Is used. Imposed constraints used In the stiff- 
ness formulation are progressively applied In a consistent mass formulation. 
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Therefore, It Is suggested that a lumped mass be used with nonconforming 
elements and a consistent mass with conforming elements. 


CLOSING REMARKS 

Progress toward reliable and efficient finite element procedures for 
rotating, flexible blade structures can best be described as uneven. A major 
reason for the slow rate of development Is the conceptual challenge facing 
blade analysts In using the finite element method to approximate the blade's 
deformation behavior under static and dynamic loads. 

To address the above challenge, blade analysts can use this manual as a 
preliminary reference tool. Comparisons are made between NASTRAN plate/shell 
and solid elements through published benchmark tests that are applicable for 
studying different blade related parameters such as twist, sweep and camber. 

The purpose of this report Is to provide a summary of blade modeling strategies 
and to document capabilities and limitations of various NASTRAN elements. The 
question of Implementing the NASTRAN capabilities for practical blades Is also 
Important In the Intended application to large displacement and frequency 
analyses, where substantially more complex finite element modeling strategies 
must be dealt with. These aspects are taken up In a companion report 
(ref. 44). 
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TABLE I. - RESULTS FOR SCORDELIS-LO SHELL 


Number of 

nnHo cnaroc 

Normalized vertical 

deflection at midpoint (B) of free edge 

HUUC o|TalC j 

per edge of 

QUAD2 

QUAD4 

QUAD8 

HEXA(8) 

HEX20 

HEX20( R) 

model 







2 

0.784 

1.376 

1 .021 

1 .320 

0.092 

1 .046 

4 

.665 

1 .050 

.984 

1 .028 

.258 

.967 

6 

.781 

1 .018 

1 .002 

1 .012 

.589 

1.003 

8 

.854 

1.008 

.997 

1 .005 

.812 

.999 

10 

.897 

1 .004 

.996 





TABLE II. - PATCH TEST RESULTS 




Maximum error In stress, percent 

QUAD2 

QUAD4 

QUADB 

HEXA(B) 

HEX20 

HEX20( R) 

Constant- stress loading 
Constant-curvature loading 

0 

30.7 

0 

0 

18 

51 .6 

0 

N/A 

0 

N/A 

0 

N/A 


TABLE III. - RESULTS FOR UNTWISTED CANTILEVER BEAM 


[Reg. = No element distortion, IRReg. = element 
distortion . ] 


Mode 

(Mesh) 

Normalized 

response at beam tip 

QUAD4 

QUADB 

HEX20 

HEXA8 

HEX20( R) 

Extension 

(Both) 

0.996 

0.999 

0.995 

0.989 

0.994 

Flapwl se 
(Reg.) 

.904 

.987 

.923 

.981 

.966 

Flapwl se 
(IRReg.) 

.752 

.971 

.967 

.746 

.927 

Edgwl se 
(Reg.) 

.986 

.991 




Edgewl se 
(IRReg.) 

.973 

.992 

.916 

.653 

.931 

Torsion 

(Both) 

.946 

.953 

1 .020 

.910 

.904 
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TABLE IV. - RESULTS FOR CANTILEVER BEAM WITH SWEEP 


Tip loading 
direction 

Normalized tip displacement In direction of load 

QUAD2 

QUA04 

QUADS 

HEXA8 

HEX20 

HEX20( R) 

In-plane (vertical) 
Out-of-plane 

m 

m 

m 

B 




R - Reduced Integration. 


TABLE V. - MORLEY'S SERIES 
SOLUTION (REF. 24) 

_3 

[w c = axlO ; 

a M = 0xlO~ 2 ; 
max 

a Mm1 n = txIO-2.] 


Angle, 

6 

Series 

solution 

a 

8 

T 

90° 




80° 

1 .4087 

4.86 

4.48 

60° 

.9318 

4.25 

3.33 

40° 

.3487 

2.81 

1 .80 

30° 

.1485 

1 .91 

1 .08 


a The principal bending 


moments are given by 


M 

M 

M 

M 


max 

min 

a 

b = 


= M + M u and 
a b 

= M - M. where 
a b 

1/2(M + M ) and 

x ^ 2 2 1/2 

tl/Mm x - M y ) 2 * M 2 y ] ,/2 


21 















TABLE VI. - PERCENTAGE ERROR 
IN Wc - BASED ON SERIES 


SOLUTION (REF. 25) 


Angle, 

6 

Mesh 

MSC/NASTRAN 

QUAD4 

TRIA3 

90° 

a 



0 

o 

00 

1 

-2.68 

.09 

.16 

-10.13 

-2.82 

-.97 

60° 

4x4 

8x8 

14x14 

-9.74 

-3.63 

-2.12 

-17.15 

-6.95 

-3.52 

40° 

4x4 

8x8 

14x14 

-24.00 

-16.26 

-12.53 

m 

30° 

4x4 

8x8 

14x14 

-29.97 
-25.93 
-21 .88 

-47.47 
-29.29 
-21 .88 
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TABLE VII. - PERCENTAGE ERROR 


IN M - BASED ON SERIES 
max 

SOLUTION (REF. 25) 


Angle, 

6 

Mesh 

MSC/NASTRAN 

QUAD4 

TRIA3 


4x4 

— 


90° 

8x8 




14x14 




4x4 

■B9 

-24.24 

00 

o 

p 

8x8 

mm 

-6.34 


14x14 

KlU 

-2.16 


4x4 

-20.54 

-25.52 

o 

o 

8x8 

-6.23 

-7.60 


14x14 

-2.54 

-3.06 


4x4 

-31.67 

-35.12 

40° 

8x8 

-15.20 

-14.59 


14x14 

-9.43 

-8.90 


4x4 

-39.53 

-42.56 

30° 

8x8 

-23.19 

-20.68 


14x14 

-13.19 

-14.50 
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TABLE VIII. - RESULTS FOR CANTILEVER BEAM WITH 

PRETWIST 3 


Tip loading 
direction 

Normalized response at beam tip 

TRIA3 

TRIA6 

HEX20 

HEXA8 

HEX20( R) 

Flapwl se 

0.993 

0.998 

0.985 

0.977 

0.995 

Edgewise 

.985 

.998 

.995 

.983 

.991 


3 Tr1angular elements used as recommended by 
NASTRAN (NO element distortions). 


TABLE IX. - RESULTS SHOWING EFFECT OF ASPECT RATIO 
(a) Aspect ratio = 1 .0 


Number of 
node spaces 3 
per edge of 
model 

Normalized lateral deflection 

at center 

QUAD4 

QUAD8 

HEXA8 

HEX20 

HEX20( R) 

2 

0.981 

0.927 

0.989 

0.023 

1 .073 

4 

1 .004 

.996 

.998 

.738 

.993 

6 

1 .003 

.999 

.999 

.967 

1 .011 

8 

1 .002 

1 .000 

1 .000 

.991 



(b) Aspect ratio = 5.0 


Number of 
node spaces 3 
per edge of 
model 

Normalized lateral deflection 

at center 

QUAD4 

QUAD8 

HEXA8 

HEX20 

HEX20( R) 

2 

1 .052 


0.955 

0.028 


4 

.991 

1 .003 

.978 

.691 


6 

.997 

1.000 

.990 

1 .066 

1 .024 

8 

.998 

1 .000 

.995 

1 .026 

1 .006 


a For elements with midside nodes, the number of 

elements per edge of model Is equal to one-half the 
number of node spaces. 
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TABLE X. - VIBRATION RESULTS FOR CAMBERED COMPRESSOR BLADE (REF. 23) 


Mode 

Symmetry 

Frequency, Hz 

Number of 

elements 

Test 
ref. 5 

01 son 
ref . 5 

2 

8 

16 

32 

1 

A 

67.5 

56.8 

86.0 


86.6/85.6 

86.6 


S 

109.6 

139.1 

136.6 

139.5 

135.5/134.5 

139.2 


S 

883.8 

279.2 

272.4 

253.7 

259 

251.5 


A 


416.6 

360.6 


351 

348.6 

5 

S 


455.3 

402.4 

413.6 

395 

393.4 

6 

A 



672.7 


564.3 

531 

533.4 

7 

A 



888.9 

818.6 

743 

752.1 

8 

S 



888.8 

867.5 

851 

746.4 


Physical properties of blade 

Sheet steel 
0.120 In. 
24 In. 

12 in. 

12 In. 

A = Anti-symmetric mode 
S = Symmetric mode 


Material : 

Thickness : 

Radius of curvature: 
Length: 

Width (developed): 






















TABLE XI. - SUMMARY OF TEST RESULTS FOR NASTRAN PLATE/SHELL ELEMENTS 
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a "Ref." = element shape not Intentionally distorted. 
a "IRReg." = element shape Is Intentionally distorted. 
a "B0TH" = Both "Reg." and "IRReg." mesh used. 
b 4 = Excellent, 3 = Good, 2 = Average, 1 = Poor, 0 = 




TABLE XII. - SUMMARY OF TEST RESULTS FOR NASTRAN 


SOLID ELEMENTS 


Benchmark test 

Mesh 3 

NASTRAN 

solid elements* 5 



— 

HEXA8 

HEX20 

HEX20(MI) 

1 . Patch test 

IRReg. 

4 

4 

4 

2. Straight beam, 

Both 

4 

4 

4 

extension 
3. Straight beam, 

Reg. 

4 

3 

3 

f lapwlse 

4. Straight beam, 

IRReg. 


3 

4 

edgewl se 

5. Straight beam, 

IRReg. 


3 

3 

f lapwl se 

6. Straight beam, 

Both 

3 

3 

3 

twist 

7. Swept beam, 

Reg. 

2 

2 

4 

edgewl se 
8. Swept beam, 

Reg. 

2 

3 

3 

f lapwl se 

9. Twisted beam, 

Reg. 

3 

3 

4 

edgewl se/ 
f lapwlse 

10. (1/4) flat plate 

Reg. 

3 

0 

1 

11. Scordells-low 

Reg. 

3 

0 

3 

shell 

12. Cambered blade 

Reg. 

3 

0 

3 

(vibrations) 
13. Twisted plate 

Reg. 

4 

3 

4 

(vibrations) 

Average r 

atlng: 

2.7 

2.4 

3.3 

Approximate grade: 

C+ 

C 

B 


Notes : 


(MI) = with modified Integration. 
a "Reg." = element shape not Intentionally distorted. 
a "IRReg. " = element shape Is Intentionally distorted. 
3,1 BOTH" = Both "Reg." and "IRReg." mesh used. 

^Rating System: 4 = Excellent; 3 = Good; 2 = Average; 

1 = Poor; 0 = N.G. 


27 












TABLE XIII. - EFFECT OF SIXTH 
D.O.F. ON DEFLECTION OF 


SHELL TIP 


Length, 

ft 

Deflection at midpoint 
of shell tip, ft 

5 DOF 

6 DOF 

Percent 

error 

25 

0.03695 

0.3696 


50 

.2374 

.2375 


100 

1 .303 

1.307 


160 

4.501 

4.509 

.18 

200 

8.064 

8.075 

.14 


TABLE XIV. - FICTITIOUS SET OF NORMAL ROTATION 
STIFFNESS COEFFICIENTS FOR TRIANGLE 
SHELL ELEMENTS (REF. 42) 



where a = arbitrary coefficient 


28 









TABLE XV. - MODAL RESULTS FOR 
SR-3C-3 PROPFAN BLADE AT ZERO 
RPM (REF. 43) 


Mode 

KR0T6 with MSC/QUAD4, Hz 


1000 

10 000 

100 000 

1 

203.5 

a ( 206 . 8 ) 
203.6 

204.1 

2 

444.8 

a ( 4 58 . 1 ) 
445.1 

446.8 

3 

663.3 

a (665.8) 

663.5 

664.7 

4 

815.2 

a (858.6) 

816.2 

821 .1 


a KR0T6 with MSC/TRIA3. 


TABLE XVI. - STRAIN GAGE RESULTS 
FOR SR-3C-3 PROPFAN BLADE 
ROTATING AT SPEED OF 
8508 RPM (REF. 43) 


Gage 

KR0T6 with MSC/QUAD4, 
In. /In. /In. 


1000 

10 000 

100 000 

1 

380.7 

a ( 402) 
380.4 

380.0 

2 

79.5 

a ( 85 . 0) 
79.4 

78.9 

3 

145.1 


146.7 


a KR0T6 With MSC/TRIA3. 
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TABLE XVII. - VIBRATIONS OF ROTATING TWISTED 
PLATES (0 DEGREES SETTING ANGLE) 

[IB = 1st bending (flap); IT = 1st torsion; 
2B = 2nd bending (flap); 3B = 3rd bending 
(flap); IE = 1st bending (edge); 


« = nondlmenslonal rotation speed.] 


Mode/n 

A 

B 

C 

D 

E 

IB 






0 

3.42 

3.4 

3.35 


3.395 

2 

4.08 

4.54 

4.4 

4.16 

4.515 

6 

7.35 

9.49 

7.34 

7.41 

9.461 

IT 






0 


21 .2 

21 .08 


20.947 

2 

21 

21.3 

21 .2 

19.56 

21.258 

6 

21 .8 

22.9 

22 

20.52 

23.51 

2B 






0 

21 .31 

21.7 

21 


21 .312 

2 

22.3 

22.4 

22 

23.03 

22.428 

6 

29.2 

29.8 

29 

29.8 

28.861 

38 






0 

59.79 

60.1 

59 


60.032 

2 

60.9 

61.3 

60.1 

62.08 

61 .164 

6 

68.9 

69.5 

68.3 

65.62 

69.468 

IE 






0 

62.23 

63.7 

63.8 


62.374 

2 

62.3 

63.8 

63.8 

62.08 

62.722 

6 

62.8 

64.3 

64.4 

64.62 

62.417 
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TABLE XVIII. - VIBRATIONS OF ROTATING TWISTEO 
PLATES (90 DEGREES SETTING ANGLE) 

[IB = 1st bending (flap); IT = 1st torsion; 
2B = 2nd bending (flap); 3B = 3rd bending 
(flap); IE = 1st bending (edge); 


n = nondlmenslonal rotation speed.] 


Mode/fi 

A 

B 

C 

D 

E 

IB 






0 

3.42 

3.4 

3.35 


3.4 

2 

4.54 

4.54 

4.49 

4.16 

4.04 

6 

9.5 

9.49 

9.48 

9.53 

9.32 

IT 






0 

20.88 

21 .2 

21.08 



2 

21.2 

21 .5 

21 .4 

19.76 


6 

23.4 

23.7 

23.6 

22.2 

21 .94 

2B 






0 

21.31 

21 .7 

20.96 


21.31 

2 

22.4 

22.4 

22.1 

23.12 

22.34 

6 

29.8 

29.8 

29.6 

30.4 

29.25 

3B 






0 

59.78 

60.1 

59 


60.03 

2 

60.9 

61 .2 

60.1 

62.15 

61 .13 

6 

69.15 

69.5 

68.5 

66.17 

69.21 

IE 






0 

62.23 

63.7 

63.8 


62.37 

2 

62.8 

63.8 

63.8 

62.8 

62.44 

6 

66.57 

64.4 

64.1 

70.64 

62.99 
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FIGURE 1. - COSMIC NASTRAN's TRIA2. 


FIGURE 2. - COSMIC NASTRAN's QUAD2 . 



FIGURE 3. - THE CLOUGH TRIANGLE FOR COSMIC/NASTRAN . 


FIGURE i|. - THE SAP SHELL 


ELEMENT FOR COSMIC/NASTRAN 




FIGURE 5. - MSC/NASTRAN's HEX8. 


FIGURE 6. - MSC/NASTRAN's HEX20. 





/ 

/ ^ 



FIGURE 9. - MSC/NASTRAN's TRIAS. 


FIGURE 10. - MSC/NASTRAN's QUAD8. 



j- SCORDEL I S-LO EXACT 
/ ANALYTICAL SHALLOW 
j SHELL SOLUTION 
/ = 3.703 IN. 


- DEEP SHELL EXACT 
SOLUTION = 3.53 in. 


O MSC NASTRAN QUAD4 
A COSMIC NASTRAN TRIA2 
□ COSMIC NASTRAN QUAD2 


FIGURE 11. - SCORDEL I S-LO ROOF. 
THICKNESS = 0.25; E = 4.32x10 8 
UNIT AREA IN -Z DIRECTION; U x 
N x N ON SHADED AREA. 


RADIUS = 25.0; LENGTH = 50.0; 
: V = 0.0; LEADING » 90.0 PER 


0 200 900 600 800 1000 1200 140 

DEGREES OF FREEDOM (1/4 MODEL) 

FIGURE 12. - CONVERGENCE OF NASTRAN ELEMENTS IN SCORDEL I S-LO 
ROOF. 


Y 



FIGURE 13. - 

t = 0.001; 


PATCH TEST FOR PLATES. 
E =1.0x10 6 ; \J = 0.25. 


a = 0.12; b = 0.24; 


FIGURE 14. - PATCH TEST FOR SOLIDS. OUTER DIMENSIONS: UNIT CUBE; 

E = I.OxlO 6 ; V = 0.25. 



FIGURE 15. - UNTWISTED CANTILEVER BEAM. LENGTH = 6.0; 
WIDTH = 0.2; DEPTH = 0.1; E = I.OxlO 7 ; V= 0.30; 
MESH =6x1; LOADING: UNIT FORCES AT FREE END. 

NOTE: ALL ELEMENTS HAVE EQUAL VOLUME. 


FIGURE 16. - CANTILEVER BEAM WITH SWEEP. INNER RADIUS = 4.12; 
OUTER RADIUS = 4.32; ARC = 90°; THICKNESS = 0.1; E = I.OxlO 7 
0 = 0.25; MESH =6x1; LOADING: UNIT FORCES AT TIP. 
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FIGURE 17. - MSGMESH TRIANGLE MESH IN MSC/NASTRAN. 



FIXED 

END 


FIGURE 18. - TWISTED BEAM. LENGTH * 12.0; WIDTH = 1.1; 
DEPTH = 0.32; TWIST = 90° (ROOT TO TIP); E = 29.0x10 6 ; 
V = 0.22; MESH = 12 X 2; LOADING: UNIT FORCES AT TIP. 


4 


1 


L 


3 


2 



t 


(A) ISOPARAMETRIC ELEMENT. 



FIGURE 19. - RECTANGULAR PLATE, a - 2.0; b = 2.0 OR 10.0; THICK- 
NESS = 0.0001 (PLATES); THICKNESS = 0.1 (SOLIDS); E = 1.7472x10 7 ; 
V = 0.3; BOUNDARIES = SIMPLY SUPPORTED OR CLAMPED; MESH = N x N 
(ON 1/4 OF PLATE); LOADING: UNIFORM PRESSURE, q * 10" 4 OR CENTRAL 
LOAD P = 4.0x10~ M . 




(C) PURE SHEAR DEFORMATION. 

FIGURE 20. - DISTORTION MODES OF STANDARD ISOPARAMETRIC ELEMENTS. 
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GAGE LOCATIONS 


TRIA3 MODEL (TRIANGLES NEARLY EQUILATERAL) 
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